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PACS 03 . 75 . Gg - Entanglement and decoherence in Bose- Einstein condensates 
PACS 05 . 45 . -a - Nonlinear dynamics and chaos 

PACS 11 . 15.Kc - General theory of fields and particles - Classical and semiclassical techniques 

Abstract. - We study the localization of bosonic atoms in an optical lattice, which interact 
in a spatially confined region. The classical theory predicts that there is no localization below 
a threshold value for the strength of interaction that is inversely proportional to the number of 
participating atoms. In a full quantum treatment, however, we find that localized states exist for 
arbitrarily weak attractive or repulsive interactions for any number (> 1) of atoms. We further 
show, using an explicit solution of the two-particle bound state and an appropriate measure 
of entanglement, that the entanglement tends to a finite value in the limit of weak interactions. 
Coupled with the non-existence of localization in an optimized quantum product state, we conclude 
that the localization exists by virtue of entanglement. 
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Spatial localization of quantum interacting particles and 
formation of bound states are of fundamental interest to 
modern physics. One intriguing aspect is the correspon- 
dence between localized states in classical and quantum 
mechanical theories [1]. Usually, one expects quantum 
fluctuations to weaken localization, as the binding of par- 
ticles with an attractive but shallow pair potential can 
be inhibited by quantum mechanical zero point motion. 
Then, localization can be interpreted essentially as a clas- 
sical property that would emerge in a quantum system 
due to decoherence [2]. On the other hand, it was recently 
suggested that localization of quantum particles may be 
achieved when they are entangled through suitable mea- 
surements [3]. Here, we consider the role of entanglement 
in the localization of specific eigenstates of a multiple bo- 
son system, e.g. the ground state. Specifically, we show 
that spatially confined interaction between atoms in an 
optical lattice induces entanglement and leads to localiza- 
tion, while the corresponding classical atomic field fails to 
localize. Remarkably, the effect that we demonstrate in 
this Letter does not depend on whether the interaction 
is attractive or repulsive. Recent experiments on forma- 
tion of repulsive atomic pairs on optical lattices [4] imply 
the possibility of experimental observation of the effect 
reported here. 



If a translationally invariant lattice with interactions is 
considered, its classical limit allows for localized solutions 
known as lattice solitons or discrete breathers [5]. A par- 
ticular realization of such a system is a Bose-Einstein con- 
densate (BEG) in an optical lattice [6]. Due to the band 
structure with Bragg reflection gaps in the optical lat- 
tice, localized soliton solutions are possible not only with 
attractive but also with repulsive interactions. Experi- 
mental evidence for the band gap solitons with repulsive 
BECs has been reported for one-dimensional lattices [7]. 
For two- and three-dimensional lattices the classical the- 
ory predicts nonzero energy and particle number thresh- 
olds for the existence of band gap solitons [8] , as opposed 
to the case of dimension one. Quantum effects in this 
system are expected to be most dramatic for a small num- 
ber of particles [9, 10]. The extreme quantum limit of a 
three-dimensional lattice has been realized in the experi- 
ment of Winkler et al. [4] , where bound pairs of repulsively 
interacting atoms have been reported using spectroscopic 
tools. According to quantum theory, these quantum soli- 
tons describe bound states of atoms that delocalize spa- 
tially [11-13]. It is an open question, whether these ob- 
served bound states persist below the above-mentioned 
classical threshold. 

In this Letter we study localization of atoms in an opti- 
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cal lattice, where interactions between atoms are present 
in a spatially confined region only. This can be achieved 
experimentally by tuning the s-wave scattering length 
by the Feshbach resonance with inhomogencous magnetic 
[14,15] or laser fields [16,17]. We show that localization oc- 
curs in the full quantum system when it is forbidden classi- 
cally. The crucial difference between the quantum and the 
classical models is the presence or absence of entanglement 
(see [18-20]) between the constituent particles. The pre- 
dicted quantum localization is due to entanglement. For 
the case of two particles we quantify the entanglement 
and show that it reaches a finite value in the limit of weak 
interactions. Beyond the specific model studied we also 
comment on the relation between entanglement in eigen- 
states, localization, and the existence of bound states for 
higher dimensional and translationally invariant systems. 
It is worth mentioning here that bound states for two elec- 
trons (fermions) in the negative hydrogen ion also appear 
only in the presence of quantum correlations beyond the 
Hartrec-Fock approximation [21]. 

The model - We initially consider the dynamics of atoms 
in a one-dimensional optical lattice in which the atoms 
interact in a spatially confined region. The Hamiltonian 
is given by 

i7 = - ^(a,\a„+i + al+ia„) + AaJaJaoOo, (1) 

n 

where ajj (a„) creates (destroys) a boson on the lattice site 
n and bosonic commutation relations [a„, aj„] = Snm hold. 
The Hamiltonian (1) describes bosonic atoms on a lattice 
that interact either repulsively (A > 0) or attractively (A < 
0) only on the single lattice site n = 0. Single atoms with 
the Hamiltonian (1) on a lattice with M sites and periodic 
boundary conditions do not localize and the eigenstates 
are plane waves l/^/M^^exp(ikn) aj,[vac) where [vac) 
is the vacuum state (no particles). However, with more 
than one particle localized states may exist around the 
site n = 0. In this current model the binding of particles 
implies spatial localization and vice versa. 

Classical treatment - The quantum Hamiltonian (1) can 
also be understood as posing a classical Hamiltonian lat- 
tice problem if we replace the particle creation and de- 
struction operators by complex valued functions of time. 
In order to enable a detailed comparison between classical 
and quantum predictions it is necessary to establish the 
precise relation between both pictures. An unambiguous 
route to relate the classical with the quantum problem can 
be found by the Hartree ansatz and variational procedure: 
For the many-body wave function with N particles we use 
the ansatz of a product state [*[f' ) = l/VNl{b'i)^\vac) 
where b'^ = l/V^X^n^n'^n creates a single particle with 
the complex amplitude Tpn on the lattice site n. The corre- 
sponding equation emerges from the standard Lagrangian 
variational procedure with g = 2X{N — 1)/N assuming 
normalized solutions with [V'nl^ — This equation 
is the discrete non-linear Schrodinger (DNLS) model with 



nonlinearity present only on the site n = 0, 
d 

iTrrtpn = -{-Ipn+l + i'n-l) + ff<5o,n I '/'O I ^'/'n ■ (2) 

ot 

This model was originally introduced to study the trans- 
port of electrons coupled to lattice phonons [22]. The 
model also applies to BECs in an optical lattice and has 
been discussed in connection with Fano resonances in the 
transport of cold atoms [23]. Here, tpn describes the com- 
plex matter wave field at the lattice site n after the intro- 
duction of appropriately rescaled dimensionless variables. 
N = l^nP is the number of atoms in the BEC. 

We have now used the Hartree procedure to derive the 
set of classical equations (2) from the quantum prob- 
lem (1). The same set of equations (2) would have also 
emerged from a more standard approach using a coherent- 
state ansatz for the many-body wave function. However, 
the Hartree procedure here serves a dual purpose in also 
characterizing the classical equations (2) as an approxi- 
mation to the quantum problem that provides strict vari- 
ational bounds for the latter. 

The model of eq. (2) supports plane wave solutions in 
the linear [g = 0) case 

tpl^"^^ = ipo exp{ikn) exp{~iLL!t) (3) 

with the dispersion relation to = — 2cosfc defining a band 
continuum [—2,2]. In addition, for non-zero g, there are 
localized solutions 

= Ae-^l"'e-'"*e''^", (4) 

with — —Ng, where the frequency [51[ ~ 2 cosh (5 > 2 
lies outside the linear band. Furthermore. 6^0 for the 
attractive interactions g < 0, where the localized solu- 
tion is the ground state, whereas 9 = n for the repul- 
sive interactions g > introduces a staggered phase pro- 
file and ipn°''^ corresponds to the highest excited state. 
From the expression = y/ N'^ — 4:/g^ for the ampli- 
tude, we find that the system exhibits a threshold for the 
existence of localized states [22], which are only found 
for N > 2/\g\. Since g may be tuned to any small 
value, the threshold for the number of particles can be 
made arbitrarily large. Conversely, for a given number of 
particles, there is a threshold value of g for localization 
to occur. Figure 1 shows the dependence of the energy 

4ai = E„ -(^'iVi'"^'* + c.c.) + f [4'°^'[4 = 1 + ^ 

on the coupling constant in the case of = 2 particles. 
In particular, no bound state is found classically in this 
system if \g\ < 1. 

Quantum case of two particles - We study the two par- 
ticle sector, where we expect to find the most obvious de- 
viations from the classical theory. In order to solve for 
the eigenstates [^f^^^) of the Hamiltonian (1) for two par- 
ticles, we introduce the projected amplitudes or two par- 
ticle wave functions (pn.m = (vac[a„am[4'(^^), which obey 
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Fig. 1: Relation between energy and coupling constant for 
two-particle defect states. The dashed line shows the classi- 
cal (Hartree) solution E^^^^ = 2(g -I- i). The shaded region 
at iJ > —4 indicates the edge of the continuum band of linear 
waves. At the classical threshold of g = 1 (the dotted line) the 
classical solution reached the continuum edge. The solid line 
shows the exact solution iS^'^' of the two-particle problem of 
eq. (8), which persist even below the classical threshold down 
to zero coupling. 

the equation 

Eipn,m = — {fn,7ii+l + <Pti,?ti-1 + V'n-l-l.m + fn-l,m) 

+ '2X6nfiS,nfi'PO,0- (5) 

This can be interpreted as the Schrodinger equation of a 
single particle on a two-dimensional lattice with a point 
defect at the lattice site (0,0). The problem is known to 
have a localized solution for any nonzero value of A [24]. 
Introducing the Fourier transform 

^k = ]j7E^^'^*''"^''"Vn,™, (6) 
m.n 

for a square lattice of M x M sites with k = (^1,^2) being 
the quasimomentum vector, eq. (5) becomes 

k' 

Looking for localized solutions with \E\ > 4 lying out- 
side the band of plane- wave energies £k = —2 (cos jjki + 
cos If ^2 ) ! we find in the limit M 00 

for the relation of the coupling parameter and the bound- 
state energy E (see fig. 1). Here, K is the complete elliptic 
integral of the first kind. It is important to emphasize 
that for 4 > > 4.05753 one has |A|, \g\ < 1 and thus no 
classical localized states persist. However, in the quantum 
case the asymptotic relation, 

E{\X\ ^ 00) ^ 2A , \E{\\\ ^ 0)1 ^ 4 + e-W , (9) 



holds and the localized state wave function is characterized 

by 

= (10) 

k ^ 

with the normahzation factor Z = -1/[RPF' (E)]. The 
bound-state wave function (/3„,m is plotted in fig. 2 for 
two classically forbidden cases. It is easy to see that this 
bound and localized state is the ground state or the highest 
energy state in the two-particle sector for attractive (A < 
0) or repulsive (A > 0) interactions, respectively. 




n 



Fig. 2: Two-particle wave function [top panels; jiySjiml, bottom 
panels: ln(|<^„m|)] for the localized state. Left panels: strong 
localization A — —0.978 and E — —4.05. Right panels: weak 
localization A = —0.779 and E — —4.01. The s-wave symme- 
try visible in the far field is a signature of entanglement, as a 
product wave function 4>n4>m would only allow four-fold sym- 
metry. Because of the (near) cylindrical symmetry of the true 
wave function, the kinetic energy cost of bringing additional 
amplitude to the site (0, 0) is significantly reduced compared 
to that of the product wave function. 

Entanglement - A system of N particles is entangled if 
the multi-particle wave function (/'ni.n2,...,njv cannot be ex- 
pressed as a product (f)\^ cj)"^^ . . . 4>^^ of single-particle wave 
functions. If the state of the system can be expressed 
by a product wave function, it is separable. The Hartree 
method becomes exact when no entanglement is present. 
Since the Hartree method is variational, it gives the best 
separable approximation in the sense that the Hartree en- 
ergy will be the closest approximation to the true eigen- 
value of the multi-particle Hamiltonian that can be ob- 
tained with a separable wave function. 

Let us discuss the two-particle problem. If the two- 
particle state were separable, due to Bose symmetry, it 
would be possible to write it in the form (t>n4'm- This is 
inconsistent with the result that in the far field, where the 
underlying lattice structure becomes less important, we 
observe cylindrical (s-wave) symmetry as seen in fig. 2. 
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A separable product approximation, on the contrary, is 
inconsistent with s-wave symmetry and is characterised 
by ridgc-hkc structures along the n = and m = co- 
ordinate axes. These structures are clearly seen in the dif- 
ference between the exact and Hartree two-particle wave 
functions shown in fig. 3. 
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Fig. 4: Entanglement in the two-particle wave function </5n,»n 
as a function of the energy E^'^^ . Above the classical threshold 
A > —1 or ij'^' > —4.05753 (shown as a dotted line) entan- 
glement is essential for localization. Shown are the von Neu- 
mann entropy S, the condensate depletion (geometric measure) 
1 — 710, and the entanglement measure V as defined in eq. 11. 
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Fig. 3: Difference vp„,m — tpn"'^'' ''Prn"^^ between the exact two- 
particle wave function and the separable (Hartree) approxima- 
tion of eqs. (4) and (10), respectively, at A = —1.26. 

In order to quantify the entanglement of the two par- 
ticles in the bound state (10) we compute various mea- 
sures of entanglement, as shown in fig. 4. First we used 
the von Neumann entropy S — Tr(plnp) [20]. Here p 
is the single-particle density matrix (SPDM) with ele- 
ments pi,j = i^^-'^(^''^)|ajai|4'(^'), normalized with F = 

^j(\I'^^^|a|ai|\l/*^^^) to have Trp ~ 1. Another measure de- 
rived from the SPDM is the condensate depletion 1 ~ uq 
(also coined geometric measure of entanglement [25]). 
Here, uq is the largest eigenvalue of p and measures the 
fraction of particles in a Bose-condensed state. Because 
p describes a pure state, I — uq measures quantum de- 
pletion, which, as we argue here, characterizes quantum 
entanglement. This would not be the case in the presence 
of incoherent, e.g. thermal, excitations. 

A third measure, V, that is amenable to analytic cal- 
culations is also shown in fig. 4. It uses projected orbitals 
defined as g„ = G~'^ Y^^ipm,n, where G = Y.,n,nVni.n- 
Since we expect for separable states that is equal to 
the product Pnffm, where gn = Vnm, the deviation 

V = ^\iPn,m- gngml"^ (11) 

is a measure of entanglement. Calculating V from eq. (10) 
analytically we find 

17 

l/-^-7H TT - 8 arctanl w 0.06858 (12) 

4 

in the limit A ^ 0, in excellent agreement with the numer- 
ical result shown in fig. 4. As this Figure shows, the differ- 



ent entanglement measures provide a similar picture, al- 
though they are in general not monotonic functions of each 
other. In particular we note that the entanglement quickly 
reaches its maximum value near the classical threshold. It 
remains finite as the two particles become infinitely weakly 
bound at A ^ 0. 

Three or more particles - We now show that bound 
states with any number of atoms larger than two exist in 
the quantum model (1) for any value of the coupling con- 
stant A 7^ as well. Without loss of generality we assume 
A < 0. We have already found a two-body bound state. 
It will suffice to show that any bound p-particle ground 
state \'ip) binds another particle for any p > 2. For this 
we have to find a (p + 1) - particle wave function [</)) with 
{(j)\H\(j)) < £;(P) -2, where E'-p^ is the energy eigenvalue of 
l^jj) and the minimum energy of a free particle is —2. We 
use the ansatz |0) = a X)n ^~'"'°'nlV')5 which is normaliz- 
ablc if a; > 1. We choose a > as a normalization constant 
to ensure = 1. We find that (01^10) < E^p'> +F{x), 

where F{x) = (1 + 2Ac)a;^ — 1 — x — x^^ + {x — x~^)p 
and c (^lajaol^) > 0. Since F(l) = -2 + 2Ac < -2 
it follows from continuity that there is an x > 1 such 
that {cj)\H\(j)) < E^P'> + F{x) < E'^p'> - 2 as required. This 
concludes the proof that bound states with any particle 
number exist in the quantum problem. 

We further remark that having found an A^-particle 
bound state that persists below the classical threshold 
A < Athresh = ^/{N ~ l)i we automatically know that 
entanglement plays an essential role in its binding. This 
is because the best separable wave function is in fact the 
Hartree approximation, which does not bind there. 

Translationally invariant systems - In interacting lat- 
tice problems with translational invariance the quantum 
eigenstates are dclocalized due to fundamental properties 
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of quantum theory. However, the existence of lattice soli- 
tons in a corresponding classical theory indicates the ex- 
istence of quantum states with local second order correla- 
tion known as quantum lattice solitons, which can be in- 
terpreted as bound states of quantum particles [9,11,12]. 
Both lattice solitons and quantum lattice solitons are char- 
acterised by frequencies and energies, respectively, outside 
of the bands of delocalized solutions in the noninteracting 
system. A framework for detailed comparison between the 
thresholds predicted by classical and quantum theory is, 
again, enabled by establishing the classical theory as a 
Hartree approximation to the quantum problem. In the 
classical theory, there is no threshold in a one dimensional 
lattice with a cubic nonlinearity (corresponding to two- 
particle interactions) but there are thresholds for higher 
dimensions [8]. The variational properties of the Hartree 
approximation guarantee that the existence of lattice soli- 
tons in the classical theory implies the existence of quan- 
tum lattice solitons but not vice versa. If quantum solitons 
exist below a classical threshold in these systems we thus 
know that entanglement between quantum particles plays 
a vital role. However, we also expect entanglement to be 
relevant for delocalised quantum soliton states above the 
thresholds (i.e. for stronger interactions). It is known that 
thresholds for quantum solitons exist in dimension higher 
than one [26] . 

Extending the current model with spatially localized in- 
teractions into more than one dimensions, there will gen- 
erally be thresholds for localization in both the quantum 
and the classical models [27]. However, these thresholds 
will generally differ. The detailed study of such systems 
lies beyond the scope of this letter and presents an inter- 
esting opportunity for future work. 

In conclusion, we have shown that localized states of 
a few atoms in an optical lattice with spatially confined 
s-wave interaction persist below the classical threshold. 
Moreover, wave function entanglement plays a crucial role 
in that localization. A one-dimensional optical lattice 
with spatially inhomogeneous interactions can be engi- 
neered with presently available techniques using magnetic 
or optically- induced Feshbach resonances [6]. Increasing 
the size of the spatial interaction domain will decrease the 
classical threshold, but it will stay finite. Thus quantum 
localization by entanglement is robust and will disappear 
only in the limit of an infinite interaction domain, where 
the classical model is known to have zero thresholds for lo- 
calized states [8] . In an experiment where interactions are 
tuned below the classical threshold the observation of lo- 
calized modes will indicate the vital role of entanglement. 
This entanglement between atoms is distillable [20] and 
could possibly be measured with entanglement witnesses 
or by reconstruction of the single-particle density matrix 
from position and momentum-space measurements. Be- 
yond the currently studied model we expect that quantum 
entanglement favors localization in other quantum lattice 
or quantum field theories as well. 
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